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Nambu-Jona-Lasinio model used to investigate low energy phenomena is nonrenormalizable, 
therefore the results depend on the regularization parameter in general. A possibility of the finite in 
four-dimensional limit and even the regularization parameter (this is dimension in the dimensional 
regularization scheme) independent analysis is shown in the leading order of the l/A^'cexpansion. 

PACS numbers: 11.30.Qc, 12.39.-x 



Ph. 

CD 



> 

o 



X 
J3 



I. INTRODUCTION 

Nambu-Jona-Lasinio (NJL) model [l| is one of the 
most popular QCD motivated effective theories used to 
understand non-perturbative low energy phenomena of 
strong interactions [l-ll] . 

Its Lagrangian contains a four-fermion interaction, an 
operator whose dimension exceeds the space-time dimen- 
sion, D, for D > 2. It is known that the four-fermion 
interaction model is renormalizable in the expan- 
sion scheme for 2 < D < 4 and the model possesses an 
ultraviolet-stable point 0-11]. It is nonrenormalizable in 
four space-time dimensions, however, the ultraviolet di- 
vergences remaining in the renormalized Green's function 
are logarithmic in the leading order of the 1/Nc expan- 
sion [3]. 

In order to construct an effective theory out of infi- 
nite number of operators involved in the NJL model one 
usually picks up operators which make a major contri- 
bution to targeted phenomena. For example, scalar type 
four-fermion operators are considered. 

As the NJL model is not renormalizable in four space- 
time dimensions some regularization methods are used 
to avoid divergences of loop integrals thereby to obtain 
finite values of predicted physical quantities. 

Three-momentum sharp cutoff regularization is widely 
used. Other regularization procedures are also stud- 
ied in NJL type models, e.g., the smooth cutoff fol. IToj. 
the dimensional regularization (DR) fll| - [l8| . the Pauli- 
Villars Il9i and the Fock-Schwingcr proper-time regular- 
ization [20l [2l| . The regularization parameter depen- 
dence is discussed in the three-momentum sharp cutoff 
and in the DR schemes in [Tsl - fisj . 

In this paper we show that for some quantities the NJL 
model behaves well in the ultraviolet limit. Moreover, 



the predicted quantities have reasonable values in the 4D 
limit; they are close to those obtained in the DR scheme 
where the dimension D < 4 is a regularisation parameter. 

The paper is organized as follows. In Sec. [ill the three- 
flavor NJL model introduced. Then we briefly review the 
dimensional regularization. In Sec. IIIIl we calculate me- 
son masses, their decay constants, etc., in the leading or- 
der of the 1/Nc expansion using DR. In Sec. lIVl the regu- 
larization parameter dependence is discussed. In Sec. |Vl 
the four dimensional limit is considered. The order of 
divergences is evaluated for n point Green functions in 
the leading order of the 1/Nc expansion. Finite values 
are obtained for physical observables. We also derive an- 
alytic relationships between observables. In Sec. IVIl we 
phenomenologically fix the model parameters and numer- 
ically evaluate the physical quantities. Some concluding 
remarks are given in Sec. IVIIl 



II. NJL MODEL 



A. NJL model 



The three-flavor NJL model including Kobayashi- 
Maskawa-'t Hooft term [i^, [13 is given as, 

'Cn.il = X] * ^ * + ^^4 + Ce, (1) 



where 



a=0 



-/f [dctgi(l-75)gj +H.C. ], 



, (2) 
(3) 



2 



the subscripts i,j are the flavor indices, i,j = w, d, s, 
and m denotes the current quark mass matrix, m = 
diag(TO„, m^, TOs). Below we consider the SU(2) isospin 
symmetric case, = m^, for simpHcity. Aq are the 
Gell-Mann matrices in the flavor space, G and K rep- 
resent the effective couphng constants for four- and six- 
fermion interaction, respectively. G and K have negative 
mass dimensions, —2 and —5 respectively, so the model 
is nonrenormalizable in four space-time dimensions. The 
determinant in £g concerns the matrix elements labeled 
by the flavor indices. We suppose the order of the cou- 
pling constants to be GNc ~ 0(1) and KN^ ~ 0(1), 
where Nc is the number of colors. 

The chiral condensates (m) generate the constituent 
quark masses, m*, inside mesons. One solves the gap 
equations to evaluate the constituent quark masses. In 
the leading order of l/N^ expansion, the gap equations 
are obtained as follows 



m* = m, + 4G(itrS'') + 2K{iivS^){iUS^), (4) 

with i ^ i ^ k ^ i. tvS'' represent the chiral condensates, 
itrS** = —{ii) which are given by the trace of the quark 
propagator. 



i trS" 



i(27r)^ 



(5) 



It is also important to discuss the mass dimension of 
the coupling constants G and K. Substituting Eq. ([7]) 
into the gap equations ^ tells us that the mass dimen- 
sions of G and K are 2 — D and 3 — 2D, respectively. 
Then the rescaled couplings are evaluated as 



G, 



M^-^G, 



2(D-4) 







K. 



(9) 



where 



It is worth mentioning that GtrS* and K{trS)^ do not 
depend on Mq, so the constituent quark masses m* are 
independent of Mq. 

This mass-rescaling parameter plays a key role to sys- 
tematically control the divergences of the loop integrals. 
We will show how the model can produce physical quan- 
tities in the ultraviolet limit. 



III. PHYSICAL QUANTITIES 

We present the prescription to calculate meson masses, 
meson decay constants and topological susceptibility in 
the leading order of the 1/Nc expansion. These quantities 
are derived from four-, two-point functions and bubble 
diagrams, respectively. 



Pion and Kaon masses 



^ — m* + ie 



(6) 



and D{= 4 — 2e) is the space-time dimensions for internal 
quark fields, "tr" in the integral denotes the trace with 
respect to the spinor and color indices. 



B. Dimensional regularization 

The quark loop integral in Eq. ([5]) is divergent in four 
space-time dimensions. One has to regularize it to ob- 
tain a finite result. The regularization dependence for 
the physical quantities are induced by this procedure. 
Using the DR we are going to take the four space-time 
dimensional limit. 

In the DR we have 



itrS' 



(2^)^/2 



r 1 



m, (m. 



*2\D/2~1 



(7) 



where the mass dimension is a function of D. The inte- 
gral like Eq. (O , needs to be multiplied by the mass scale 
parameter A/q ~^ in order to correct the mass dimension. 
Note that in the previous studies [l6l - [T8| , this parameter 
is called "the renormalization scale". However, in this 
paper, we name it "mass scale parameter" to clearly dis- 
tinguish our treatment from the renormalization. Then, 
the "rescaled" chiral condensates are given by 



M, 



n). 



(8) 



The masses of pion and kaon arc obtained by evaluat- 
ing the poles of their propagators 



1 - 2KcUp{k^)' 



(10) 



using the random-phase approximation and the l/Nc ex- 
pansion, where a labels the channel isospin and P de- 
notes the meson species. The explicit form of the flavor- 
dependent effective couplings is given by 



K3 = G+^K{itrS'), for 7r°, 



(11) 



KeEEG+ -K{itYS'^), for K", K". (12) 



Hp is the meson self-energy. 



n,(fc2) = 2nr(fc'), 
nK(fc2) = 2nr(fc2), 



(13) 
(14) 



where IC^ {k ) is the following loop integral: 
K'ik')^ J -^tr[j,S\p + k/2)^,S^p-k/2)] 
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trl 



' I i{2Tr)D (p2-mp)[{p-k)'^-m*■ 



il6) 



The trace runs over color and spinor indices, then trl = 

The conditions which determine the pion and kaon 
masses are 



l-2KenK{ml)=0. 



(17) 
(18) 



These equations are the relations between model param- 
eters and the input physical quantities and tok- 



B. rj and ?/ masses 

rj and rj' mesons are the mass cigcnstatcs for rj^ and 
TjQ mixing. In the random-phase approximation and the 
1/Nc expansion, the propagator of the r/ — ?/ system is 
given by [3|, i^] 



2K^ 



i^2K+n{k^)' 



(19) 



where the effective coupling and the self-energy 11 
are the 2x2 matrices 



K+ = 

n =^ 



Kso Kss 



(20) 
(21) 



with 



A'oo = G - -/V (i trS*" + 2i trfi-"), 

O 

6 

Ko8 = Kso = -^/v(^tr5^ -itiS"), 

D 



and 



noo(fc2) = -[2nr(fc2) + nr(fc^)]. 



no8(fc2) = n8o(fc' 



2^2 



[nr(fc')-nr(fc2)] 



The masses of 77 and rj' are obtained by solving 

det[l - 2i<:+n(m|)] = 0, (22) 

where mp denotes the mass of i] or r/. They are also 
obtained via diagonalization of Eq. ([T^ 0, [25| . 



C. Pion and Kaon decay constants 

The decay constants of pion and kaon, /p are defined 
by the matrix element of axial current between the meson 
and vacuum states, 



-M, 



4-D 



(27r)^ 



tr 



T 

7m75^5"5p,,75T]S'^ 



, (23) 



where the meson-quark-quark effective coupling gpqq is 
defined by 



(3: 



Pqq) 



M, 



4_^9np(fc2) 



(24) 



For notational simplicity, when it is obvious, we will omit 
the subscript "rs" indicating that the quantity is rescaled 
with the help of a power of Mq. In the leading order 
of 1/iVc expansion, the decay constants are calculated 

as m. 



M, 



4-D 



X tr / dx 



KlusiO) + {ml - ml) 
d^p X 



,0 J i(27r)D {p2^Lus{0)+ieV 
where Jus is defined by 

Jns(fc') = (fc') + 2{{m* - m:)2 - 

"1 '■ d^p x{l - x) 



(25) 



(26) 



X tr / dx 



with 



Lij[k'^) 



i{2^)D {p^~Lus{k-^) + ieY' 



{mf - mf)x - k'^xil - x). 



(27) 



Equation (pSj) is used to fix the mass scale parameter Mq. 



D. Topological susceptibility 

The topological susceptibility % is defined by the corre- 
lation function between the topological charge densities, 
Q{x), at different points 

X = y d^x(0|TQ(x)Q(0)|0) 

connected.! 

(28) 

where 



Qix) 



32n 



.papa^.. ^ 2/i Im[det q{l - 75)9] , (29) 



g is the strong coupling constant of QCD and is 
the field strength for gluons. Equation ([^5]) should be 
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multiplied by Mq ^ to adjust the mass dimensions. In 
of] 

(itr5")2 



the leading order of l/iVc expansion x is given by [2; 



X 



M, 



D~4 ' 



(itrS"')(«trS"') \ 



1 



-=(2ztr5^ +ztr5")(noo(0),no8(0)) 
+ (ztr5^ - ztr5") (Hog (0) , Uss (0)) | A+ (0) 

x{-L(2.tr5. + .«-,(n«»<°> 



-L,,..s--..5-)(n»(°) 



(30) 



Thus X is evaluated by the quantities already obtained 
above. 



IV. STRATEGY 

As mentioned in the introduction, our goal is to test 
the regularization parameter dependence (or indepen- 
dence) in the NJL model. 

We use the abstract denotations: m(e M) are the 
model predictions, p(G V) are the model parameters, and 
i(S I) are the input physical quantities. The capitals 
M^V and X represent the sets of these quantities. The 
model predictions can be performed if all the parameters 
are known, so the model determines some function T^p, 



m 



(31) 



where the superscript R indicates a regularization proce- 
dure. In the similar manner one can obtain the functions 
J-^ which connect V and I through the parameter fit- 
ting. Thus, the model relates I to Al, as 



m«(I)=J-«i(I). 



(32) 



The resulting values of m should not depend on regular- 
ization methods, if they correctly capture the physics in 
question. In the next section, we obtain these functions 
J-^i, J-«pand 

Our model with the DR has six free parameters 
= {wu, TTis, G, JiT, D, Mo}. Four of them, V4 ~ 
{msi G, K, Mq}, may be fitted to the input meson prop- 
erties, X4 = {m7r, /tt, tok, ^r;'}- After this partial 
parameter setting, predicted meson properties, A4 ~ 
{/k7 "T-r;, X: ■ • ■ } Written as the functions of the re- 
maining parameters niu and D, 



m 



(33) 



We cannot obtain finite values for all predicted quantities 
in the 4D limit, all of m are finite at D < 4. However, it 
is possible to obtain finite values for some of m at the four 
dimensional limit. Furthermore, some of physical quan- 
tities may be independent of the parameter D [l^ . We 



are interested in this aspect of the NJL model which does 
not depend on the regularization parameter. Namely, it 
is interesting that in the model discussed here there are 
theoretical relations between some observables which are 
finite (although the model is not renormalizable) and, 
moreover, they do not depend on the regularization pa- 
rameter. We mean the relations between input and pre- 
dicted observables derived with the help of the fitting 
parameters procedure. 

Once we put the actual numbers into I4 and fix D, 
m become the functions of niu- In particular, we shall 
evaluate them in the — >■ 4 limit, then compare between 
each other results obtained in the three regularization 
ways: (1) DR with D 4, (2) DR with fixed D by 
X, and (3) three-momentum (3M) cutoff method. These 
regularization schemes can symbolically be written as 

(1) m4D(m„) = T^^{m.a)\v,^i,M^i, 

(2) mi^«(TO„) = J-™(mO|{P„i5}^{i„x}, 

^3U/ \ _ T-3M/ 



where Vf^ = {nis,G,K,A}, with the momentum cutoff 
A. Note that the total number of the parameters in the 
cutoff method is five [l^ : since this method does not need 
the mass parameter. 



V. FOUR DIMENSIONAL LIMIT 

We demonstrate the calculation in the four dimen- 
sional (4D) limit and express the predicted physical quan- 
tities in terms of the input meson properties. By virtue 
of taking the 4D limit, the relations among the physical 
quantities can be simplified through the leading order of 
e expansion. It is even possible to obtain analytic expres- 
sions for some quantities, which will be discussed below. 

A. The summary of the Sec. IIIII 

To make the model prediction we need to know the 
following six quantities: 

m*, m*, G, K, m^, Mq. (34) 

These are evaluated by the gap equations 

ml = niu + AGiitrS") + 2K{itrS''){itiS'), (35) 
=m,, +4G(itr5")+2i^(itr5")2, (36) 

and the conditions for meson properties 

l-2A-3n,(77i2) = 0, (37) 

l-2KeTlKiml) = 0, (38) 

det[l - 2K+Il{ml,)] = 0, (39) 



(40) 



In the following, we will calculate the quantities by 
solving these six equations in the 4D limit. 
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B. Constituent quark masses m*, m* 

The constituent quark mass to* is obtained from 
Eqs. ([55]) and (P7|) . The gap equation ([55]) can be rewrit- 
ten as 



wJl = m„ + 4/s:3(itrS'"), 



(41) 



which in combination with Eq. ([57)1 enables one to get 
rid of -ftTs and thereby to calculate m* without knowing 
the values of G and K. 

In the 4D(e — >■ 0) limit, the chiral condensate itrS*' and 
the self-energy loop integral IV^ {k'^) can be expanded in 
powers of e as 



itrS" 



N, 



(42) 



n5'(fc') = ^ [k' ~ 2(TOf + TOf - m*m*)] . (43) 
After some algebra we arrive at 



(44) 



Here we adopted the negative solution of m*, since the 
positive solution is unstable. 

We thus find the analytic expression for to* as a func- 
tion of TO,r- On the other hand we apply numerical anal- 
ysis to obtain to* by simultaneously solving Eqs. (j37| . 
and 




Thus the couplings are the functions of to*, G{m*) and 
Kijnl). With the help of Eq. ([39]), we can determine to* 
and substitute it in G{ml) and K{m*) to get G and K. 

Note that G and K are of the order of e and e^, respec- 
tively. In the naive e — limit G and K tend to [24| . 
However, the combinations GtrS* and i4r(trS')^ appears 



to be finite because the chiral condensate, Eq. (|42p . is of 
the order of e~^. Then, the dynamical masses, Eq. (U), 
can have larger values than even in the limit of e — > 0. 

From the order counting of 1/e for Eqs. (|35|) . ([36]) and 
(P^. we find the following relation. 



L-Ng-2Nk 



l-n/2 



(49) 



where L, Nq and A'^x arc the number of quark loops (chi- 
ral condensates), G and K, respectively, n is the number 
of the external quark lines (n = 2^, / = 0, 1, 2, • • • , Z are 
the number of external meson lines) . Then n point func- 
tions become finite except the bubble diagrams. 

For instance, the above discussion is manifested in the 
expressions for the pion and kaon propagators in the limit 
D ^ 4, 



47r^e 



1 



(50) 



These are derived through the substitution of the ob- 
tained G, K and Hp into Eq. (fTO)) . 



D. Current strange quark mass 



C. Couplings G, K 

The couplings G and K can be written as the functions 
of TO*, TO*, TOtt and tok. 

From Eqs. ([57| and one derives the following re- 
lations 



G 
K 



-itrS'"nK + itrS-^n^ 
2(itrS'* - itrS'")H^HK ' 

Hk-H^ 
{itrS' - itrS'«)H^HK ' 



(45) 
(46) 



where we use the abbreviated notations of 11-^ ~ H7r(TO^) 
and Hk = Hk(to|). By using Eqs. and (gS]) 

we arrive at the following expressions for the couplings 



G 



K 



2^2e 



1 



m% — 2to*2 
leTT^e^ 1 



TO^ — 2(771*^ + TO*^ — m^m*) 
0(e2), (47) 
1 

TO.|- — 2(to*-^ + TO*^ — ml^m*) 



1 



— 2to*- 



0(e=^ 



(48) 



We have seen that to*, to*, G and ii' are determined 
from Eqs. dSS]), dSZ]), ® and Substituting these 

into Eq. p6p . it is easy to evaluate rris. The numerical 
result will be shown in Sec. IVIBI 



E. Mass scale parameter Mq 

In this subsection, we shall set the remaining parame- 
ter Mq, the mass scale, by using the pion decay constant 
in Eq. go]). 

In the 4D limit, the integral /„„ (Eq. HI])) becomes 
lim/„„==hm^^-^r(e)y^ dxL,- 



lim 



Hence we obtain the mass scale from Eq. ([?5]): 

4^26 



lim Mn = lim 



£->0 



>0 iVr TO 



*2 ■ 



(51) 



(52) 



Since to* is the function of to^t, (JH]), the value of Mq is 
determined by to^ and /^. Mq goes to in the limit of 
0. 
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It is worth mentioning that by using Eqs. ((52)) . ([26 
and (P7)) . /k can be written as 



f2 _ r2 
■Ik — Jit 



4m*2 



(53) 



since Jij coincides with /y in the leading order of the e 
expansion. 



F. Chiral condensates {uu), (ss) 

From the Eq. (^^ the chiral condensates (m) are of the 
order of e^^, so they diverge in the e — > limit. However, 
the rescaled chiral condensates {ii)rs (= M^~^ [ii)) are 
finite, since the order of Mq~^ is e as seen in Eq. ([5^ . 

The form of the rescaled chiral condensate ('Sw)is is 
obtained from Eqs. (|42|), dM]) and ([52]) as 



(Mlt)i.s = - 




(54) 



If we ignore the last term of Eq. (|54p . the above equa- 
tion coincides with the Gell-Mann-Oakes-Renner rela- 
tion [2i, [2i|. With the help of Eq. one has the 
analytic expression for {ss)^s- 




AjTli, 



In the limit Jk ^ fir, {ss)i-s coincides with 



(55) 



VI. NUMERICAL RESULTS 

To evaluate the physical quantities discussed in Sec. lIIII 
we employ the following input meson properties: 

= 138MeV, U = 92MeV, 
rriK = 495MeV, m^, = 958MeV. 

We first evaluate the constituent quark masses and the 
strange current quark mass, then calculate the meson 
properties /k, m^, and the topological susceptibility x- 

As mentioned in Sec. II VI we compare the results ob- 
tained in the the three cases: (1) DR with 4D, (2) DR 
with D{x), and (3) 3M cutoff method. In the second 
case, we select the value x^^* = ITOMeV [11] and call it 
simply the DR results. 



A. Constituent quark masses m*, m* 

Before evaluating the observed quantities, we consider 
an unobserved quantity, the constituent quark mass. In 
Fig. [U TO* and m* are shown as the functions of to„ in 
the range 3MeV < m„ < 6MeV. Since to„ is contained 
in the denominator of Eq. (j44p , to* depends strongly on 



-1500 F 



> -2000 



-2500 
-3000 



3.0 3.5 4.0 4.5 5.0 5.5 6.0 
m„ (MeV) 

FIG. 1. and m* as the functions of m„. The solid and 
dashed lines are the results of m* and m* , respectively. 



TO„. Note that to* and to* can be positive according 
the discussion in Ref . . The values of |to*| and |to*| 
are considerably larger than in the frequently used 3M 
cutoff case, m* - 300McV and to* - 500MeV. There- 
fore a large regularization dependence is found for an 
unobserved quantity. Note that to* can be obtained by 
solving Eq. ([53]) analytically. However, this solution does 
not lead a realistic value of m^/ . In other words the ob- 
served value of the kaon decay constant, /k, can not be 
consistent with the realistic value for to,,' . 



B. Current quark mass nis 

The plots of TUs obtained in the 4D, DR and cutoff 
cases are shown in Fig. [2] We see that TOg increases 




4.0 4.5 5.0 
m„ (MeV) 



6.0 



FIG. 2. The strange quark mass nia as the function of rriu. 
The black line shows rris in the 4D case. The circles and 
squares are DR and cutoff regularization results, respectively. 
The gray bound is the experimental region. 

linearly with respect to to„. It is interesting to note that 
the DR and cutoff results show behavior similar to the 
4D result. The results around to^ ~ 5MeV cross the 
experimental region which is evaluated at IGeV [l^. 
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C. Kaon decay constant /k 

Figure ini displays the results of /k with the experimen- 
tal value = llOMeV. The resulting /k in the 4D case 



115F 



110 



^ 105; 



100 : 
95 : 

90 L_ 
3.0 



3.5 4.0 4.5 5.0 

m„ (MeV) 



5.5 6.0 



FIG. 3. The kaon decay constant /k as the function of m^. 
Black and gray lines are the obtained /k and the experimental 
value. The circles and squares are DR and cutoff regulariza- 
tion results, respectively. 

is smaller than its experimental value. The DR plots are 
similar to the 4D case, and they are a few MeV closer to 
the experimental line. On the other hand, /k decreases 
with increasing m„ in the cutoff rcgularization, which is 
the opposite to the 4D and DR cases tendency. The cut- 
off results for small niu region receive the large effect of 
the parameter dependence. 



D. Eta meson mass m„ 



We put the result of in Fig. |4l In the 4D case, 



> 



560 1 
540 
520 
500 



480 E* 

460 
440 
420 



3.0 3.5 



4.0 4.5 5.0 
m„ (MeV) 



5.5 6.0 



FIG. 4. The eta meson mass m,, as the function of m„. The 
black and gray lines are the obtained rrir, and the experimental 
value, m™ = 548MeV. The circles and squares are the DR and 
the cutoff results. 

is around SOOMeV at = S.OMeV, and it slightly in- 
creases with respect to TOk. For all the region, the values 
are smaller than in the experimental data. Contrary to 
the results seen in /k, both the DR and cutoff cases are 



worse than the 4D case in terms of comparison with the 
experimental data. 



E. Topological susceptibility x 
In Fig. O X is calculated as a function of . One sees 
190 r 




4.0 4.5 5.0 

OT„ (MeV) 



FIG. 5. The topological susceptibility x a-s the function of 
niu- The black line is the obtained x in the 4D. The circles 
and squares are the DR and cutoff results, respectively. The 
gray bound is the result of lattice simulation [301 ]. 

that, in the 4D case, x becomes larger when to„ increases. 
The DR result is trivially fixed at x^^^ = 170MeV, be- 
cause X is the fifth input quantity in this case. The 
gray bound shows the results of lattice simulation jsoj . 
X^/^ = 170±_r 174±7MeV, and Witten-Veneziano mass 
formula 0,111, x^^^ = 179MeV. The 4D result is close 
to the values of lattice simulation. The 4D and DR re- 
sults arc plotted inside the lattice region. However, al- 
most all squares are located outside of the lattice region, 
so the cutoff case is worse than the 4D and DR in this 
context. 



F. Chiral condensates (uu), (ss) 

The chiral condensates are shown as functions of m„ 
in Fig. [6] where we use the obtained /k to evaluate 
Eq. ([55]) . The rriu dependence can easily be read off 
the explicit forms in Eqs. ([54]) and (|55p . It is interest- 
ing to note that, although the values of |m*| in the 4D 
and DR methods arc much larger than in the cutoff case, 
TO* ^ ©(lOOMeV), the chiral condensates almost coin- 
cide in the 4D and DR regularizations. 



VII. CONCLUSION 

We have constructed a NJL model with DR which is 
convergent in the 4D limit, and found that it is possible 
to control the ultraviolet divergences within the frame- 
work mentioned in the section |Vl The obtained results 
describe meson properties pretty well. We believe that 
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FIG. 6. The chiral condensates {uu)rB and (ss)rs as the func- 
tion of rriu- The sohd and dashed lines indicate (uu) rs and 
{ss)rB, respectively. 



the treatment prescribed here has possibilities for future 
work. 

(I) We have shown that the model is free of divergenees 
by virtue of introducing the mass scale parameter. The 
mass scale runs according the dimension in the dimen- 
sional regularization scheme where the ultraviolet limit 
corresponds to 4D limit. We find that it is possible to 
remove all the divergence within the model frame work. 
Here the mass-scaling determines the scale of the model, 
so this is reminiscent of the renormalization flow. We 
have explicitly shown the ultraviolet behavior via numer- 
ical calculations. 

(II) It is a non-trivial question whether so constructed 
model can produce reliable results, because important 
contributions may be dropped due to taking the 4D limit. 
However the obtained results indicate that the model be- 
haves well even in this limit. 

(HI) The meson properties are described without the 
regularization parameter. The model predictions do not 
depend on D due to the effect of the mass rescaling. This 
indicates that we can consider the ultraviolet limit even 
in the nonrenormalizable model. 

(IV) We found that, in this model framework, it is pos- 
sible to obtain analytical relations between meson prop- 
erties and the chiral condensates. This point has practi- 
cal importance as the model calculations are significantly 
simplified. In particular, Eq. (|50[) has quite simple form 
for the meson propagators. This form implies that we 
can introduce meson propagators in the same form as it 
appears in the usual perturbative expansion techniques. 
Therefore the method is expected to be useful in more 
complicated problems such as for example in the three- 
body formalism. 

We have also clarified the following less important 
points: 

(V) By applying the effects of the mass scaling, we have 



clarified the role of the mass dimension. It is known that 
the couplings should become zero in the ultraviolet limit 
so that the theory well behaves in the ultraviolet limit . 
We have effectively incorporated the "running" couplings 
by using the mass rescaling to control the divergences 
with keeping the finite values of the rescaled couplings. 

(VI) Related to the point (I) we have introduced the 
effective coupling scale in the couplings G and K, and 
this works successfully. This is trivial if one deals with 
the renormalized theory where the ultraviolet behavior 
is well known through the renormalization group flow. 
However, in the NJL model, the background gluon contri- 
butions are implicitly expressed in terms of the effective 
coupling strength. Therefore the couplings regularized in 
specific ways contain important dynamical information of 
the model. Our results indicate that the model with the 
DR does not miss dominant physical contributions when 
we take the high-energy limit. 

The model predictions are intimately related to the 
employed regularization prescriptions, because they in- 
clude the background dynamical information as was ex- 
plained in the point (VI). Then, it may be interesting to 
investigate whether the ultraviolet behavior of the model 
different regularization schemes leads to the similar con- 
sequences. We believe the DR method is especially good 
in this context, because it is expected to preserve the 
required symmetry of the model just as it does quite suc- 
cessfully in the formal perturbative quantum field theo- 
ries. This statement is confirmed by phenomenological 
results of the paper. 

Concerning the point (II) we also found that the phys- 
ical quantities calculated in the DR with D < A and in 
the cutoff regularization have values similar to those ob- 
tained in the DR at the 4D limit. 

However, there is a discrepancy between the model re- 
sults and the experimentally observed values. We believe 
that this discrepancy should be explained by introducing 
higher dimensional operators. 

Calculations are drastically simplified if one uses our 
method thanks to taking the 4D limit. Thus we hope 
that the regularization parameter independent approach 
discussed in this paper can be useful to systematically 
introduce higher dimensional operators and to calculate 
higher order corrections. 
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